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MASS AMBIGUITIES IN CASCADE DECAYS 
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We review the use of invariant mass distributions in cascade decays to measure the masses of New 
Physics (NP) particles in scenarios where the final NP cascade particle is invisible. We extend 
earlier work by exploring further the problem of multiple solutions for the masses. 
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1. Introduction 

Inspired by the precise determination of 
the dark matter relic density made possible 
with data from the WMAP satellite 1 ' 2 - 3 New 
Physics (NP) models with a weakly interact- 
ing massive particle have flourished over the 
last years. Such particles, whether they ap- 
pear in models of supersymmetry, universal 
extra dimensions or little Higgs theories, are 
invisible to the detectors of collider experi- 
ments, and to be viable dark matter candi- 
dates, certainly stable on the scale of the cur- 
rent age of the universe. The conservation of 
a NP quantum number, be it R-parity, KK- 
parity or T-parity, ensures a cascade decay of 
all heavier NP particles into the lightest NP 
particle. Measuring NP particle masses in 
such scenarios by reconstructing mass peaks 
is then in general no longer possible. 

The well studied alternative to this is to 
use measurements of the endpoints of invari- 
ant mass distributions. 4,5 ' 6,7 ' 8,9 ' 10 As we can 
express endpoint positions in terms of the 
NP masses, we can in principle solve for the 



masses if we have at least as many endpoints 
as unknown masses. However, there are non- 
trivial problems with this method. Generally 
the expressions for the endpoints are com- 
plicated functions, leading to possible mul- 
tiple solutions, and strong correlations be- 
tween the measured masses. 11 As an illustra- 
tion of this problem, we consider the decay 
chain 



and show in Fig. 1 the quantity 
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for a scan over the SPSla 12 mSUGRA scalar 
(too) and gaugino (mi/2) mass plane for a 
trilinear coupling Aq = —too, tan/3 = 10 and 
/i > 0. The sum in (2) runs over the squark, 
slepton and neutralino masses in the decay 
chain (1). 

We find that for virtually all points in 
the mass plane where the decay chain (1) ex- 
ists, and where the lightest supersymmetric 
particle (LSP) is not the stau (at low mo), we 
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Fig. 1. Scan of fi, as given by Eq. (2), over the 
SPSla mass plane. In the top white area the de- 
cay chain (1) is not kinematically accessible, while 
the bottom white area has a charged stau LSP. 

get a false solution from the set of endpoints. 
Along the middle of the physical wedge in 
this plane, and intersecting the SPSla line, 
the above quantity can become rather large. 
In an experiment one might also encounter 
the problem of "feet" , i.e. structures in the 
distributions near the endpoints that are hid- 
den by background. 8,13 

We have previously started an investiga- 
tion into the possible improvements to this 
method that can result from fitting the whole 
invariant mass distribution to an analytical 
expression derived for its shape. 13 We found 
that knowledge of the shape will predict any 
possible irregular features at the endpoints 
for a measured set of masses, thus essentially 
removing the problem of feet. Here we want 
to investigate further what effect knowledge 
of the shape of the invariant mass distribu- 
tion has on the problem of multiple solutions. 

2. Multiple solutions 

Using as our example the four possible a in- 
variant mass distributions for the SUSY de- 

a We can form the invariant mass of the two leptons, 
the two leptons and a jet, but assuming that we can- 
not, in general, identify which lepton is the lepton 
nearest to the squark in the decay chain, we use jet- 
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Fig. 2. Invariant mass distributions for the SU1 
benchmark point and the decay chain (1) with ar- 
bitrary normalisation. In blue (solid) are the distri- 
butions for the nominal masses, in red (dashed) are 
the distributions for the masses in the false solution. 

cay chain (1) we show in Fig. 2 the shape of 
these distributions, as found in Ref. 13, for 
the nominal masses of the ATLAS SU1 stau- 
coannihilation benchmark point. (One may 
construct other invariant mass distribution 
by imposing cuts, but for simplicity we will 
neglect them here.) 

For SU1 there is also another set of 
masses that give the same endpoints for the 
invariant mass distributions. If only end- 
points are measured we have two indistin- 
guishable solutions for the SUSY masses. 
The two sets are given in Table 1 as the nom- 
inal masses of the benchmark point, and the 
false solution from the nominal endpoints of 
the benchmark point. 

In Fig. 2 we see that there are only small 
visible differences in the shapes between the 
two mass sets of SU1 for three of the distri- 
butions, but for m g ;(i ow ) there is a marked 



lepton combinations that give the highest and lowest 
invariant mass. 
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Table 1. Nominal and false sets of 
masses for benchmark points [GeV]. 



Set 


m-o 


rar 


m-o 


m 9L 


SU1 nom. 


137 


254 


264 


760 


SU1 false 


122 


127 


246 


744 


SU3 nom. 


118 


155 


219 


631 


SU3 false 


347 


411 


452 


900 



change. Similar differences in the shapes 
can also be seen in three of the distributions 
for SU3, the bulk-region benchmark point, 
shown in Fig. 3, although they are individu- 
ally somewhat less distinct than for SU1. For 
both cases it seems likely that the false so- 
lution could be ruled out, with the expected 
high statistics at the LHC for SUSY scenar- 
ios with sub-TeV squark masses. For the 
SU3 benchmark point the large differences in 
mass between the solutions should also show 
up clearly in the measured cross sections. 

This is unfortunately not the situation 
for all models. For the SPSla bulk-region 
benchmark point 12 the differences in the 
shapes are very small. After adding smear- 
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Fig. 3. Invariant mass distributions for the SU3 
benchmark point. See caption of Fig. 2 for details. 



3 

ing from detector effects, combinatorics from 
picking the correct jet, initial and final state 
radiation etc., separating the two mass sets 
for SPSla will be very difficult and require 
very large statistics and an extremely good 
understanding of the detector, if at all possi- 
ble. 

It is then natural to ask how effective 
the shape can be in distinguishing multiple 
solutions over a wide range of SUSY models. 
In Fig. 4 we show a scan over the SPSla mass 
plane. We plot the value of 

1 /* m i IlaX 

C = ^E \fi(m)- 9i (m)\dm (3) 

where fi and gi are the invariant mass distri- 
butions for the false and nominal sets of mass 
solutions (normalized to J dmf(m) = 1, 
J dmg(m) = 1), with the i-index running 
over the number of distributions N (here 
N = 4), and where mf ax is the endpoint 
of the i-th distribution. The numerical pre- 
factor ensures that D < 1. 

The difference in shape between the 
nominal and false solutions, as quantified 
by the D-value is fairly constant at around 
D « 0.07, except for a wedge in the mass 
plane, where we have a transition between 
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Fig. 4. Scan of D, as given by Eq. (3), over the 
SPSla mass plane. In the top white area the decay 
chain (1) is not kinematically accessible, while the 
bottom white area has a charged stau LSP. 
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different forms for the analytical expressions 
for the m q i(i ow ) and m 9 /(hi g h) distributions. 13 
In this region the D-value changes rapidly 
from low to high values, and from Fig. 4 
we see that the SPSla benchmark point suf- 
fers from an unfortunate position, giving it 
a value of D = 0.016. Comparing to Fig. 1 
we find that while large mass differences can 
give large differences in shape, there are cer- 
tainly areas where the correct solution could 
be readily identified from the shape, even 
when the mass differences are small, and thus 
where a comparison of cross sections would 
be difficult. 

Due to smearing from the experimental 
effects mentioned above, it is difficult to cal- 
culate exactly for what D-value two sets of 
solutions can be distinguished for some given 
integrated luminosity at the LHC. However, 
assuming that the two distributions of SU1 
are distinguishable, and using as a basis for 
comparison that we find a D-valuc of D = 
0.070 for SU1, this suggests that we should 
be able to distinguish between the nominal 
and false solutions in most of the parameter 
space shown in Fig. 4. We have performed 
such scans for different values of Aq and 
tan (5, in the mass planes of the SPSlb, SPS3 
and SPS5 benchmark points, 12 and find sim- 
ilar results there. 

3. Conclusions 

We have discussed aspects of the use of an- 
alytical expressions for invariant mass distri- 
butions in determining masses in NP scenar- 
ios, using a well studied SUSY decay chain 
as our example. We have focused on the pos- 
sibility of using differences in the shapes of 
distributions to remove false solutions that 
appear if only the endpoints of the invari- 
ant mass distributions are used to determine 
the NP masses. Our results from scans of 
the mSUGRA parameter space indicate that 
with the exception of certain limited areas, 
the shapes of the distributions should in gen- 
eral be well suited to reject the false solu- 



tion, given that enough statistics are avail- 
able, and may be applied even when a com- 
parison of cross sections would be difficult. 
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